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Atmospheric neutrino and solar neutrino data from the first phase of Super-Kamiokande 
(SK-I) are presented. The observed data are used to study atmospheric and solar neutrino 
oscillations. Zenith angle distributions from various atmospheric neutrino data samples are 
used to estimate the neutrino oscillation parameter region. In addition, a new result of the 
L/E measurement is presented. A dip in the L/E distribution was observed in the data, as 
predicted from the sinusoidal flavor transition probability of neutrino oscillation. The energy 
spectrum and the time variation such as day/night and seasonal differences of solar neutrino 
flux are measured in Super-Kamiokande. The neutrino oscillation parameters are strongly 
constrained from those measurements. 



1 Introduction 



Super-Kamiokande is a 50,000 ton water Cherenkov detector located 1,000 m (2,700 m water 
equivalent) under Mt. Ikenoyama at Kamioka Observatory, Gifu Prefecture, Japan. The de- 
tector is a cylindrical tank and is optically divided into two regions. The inner detector (ID) 
is instrumented with 11,146 inward facing 20 inch PMTs which give a photo cathode coverage 
of 40 %. The outer detector (OD) completely surrounds the ID with the thickness of 2.05 m to 
2.2 m water and is monitored by 1,885 outward-facing 8 inch PMTs. The OD works as a veto 
counter against cosmic ray muons. The charge information observed in the OD is also used to 
separate event sample in the L/E analysis. 



2 Atmospheric neutrinos 

2. 1 Introduction 

Atmospheric neutrinos were observed in Super-Kamiokande during a 1489 hve day exposure 
which corresponds to 92kiloton-yr. The atmospheric neutrino events are classified into fuUy 
contained (FC), partiahy contained (PC) and upward going muons. The vertices of neutrino 
interactions are required to be inside the fiducial volume of the ID for FC and PC events. If 
the tracks of entire particles are contained inside the ID, the event is classified into FC. While 
one of the particles, mostly muon, exits the ID and deposits visible energy in the OD for PC 
events. Each observed Cherenkov ring is identified as either e-like or //-like based on the ring 
pattern. The directions and the momentum of charged particles can be reconstructed from the 
ring image. Upward going muons are also observed as an atmospheric neutrino sample. These 
muons originate from high energy neutrino interactions with the rock surrounding the detector. 
Upward going muon events are classified into upward stopping muons having only an entrance 
signal in the OD, and upward through-going muons having both entrance and exit signals. The 
atmospheric neutrino events in Super-Kamiokande are predicted by a detailed Monte Carlo 
simulation—. 

2.2 Zenith angle analysis 
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Figure 1: Zenith angle distribution for fully-contained single-ring e-like and /i-like events, multi-ring /i-like events, 
partially contained events and upward-going muons. The points show the data and the solid lines show the Monte 
Carlo events without neutrino oscillation. The dashed lines show the best-fit expectations for i^^ <-^ Vr oscillations. 



Figure^shows the zenith angle distributions of each atmospheric neutrino sample. The data 
are compared with the MC expectation. The observed muon neutrino events exhibit a strong 
zenith angle dependent deficit compared with the expectation without neutrino oscillation. On 
the other hand, the observed electron neutrino events are consistent with the prediction within 
about 10 % normalization. 



A fit to the FC, PC and upward-going muon data is carried out assuming 2-fiavor f„ 



<-5- Vt 



oscillation. In 2-flavor v.. 



P{^, 



v-T oscillations, survival probability of v^ is expressed by : 
'l.27Am2(eV2)L(km)\ 



1 - sin^ 261 sin^ 



^(GeV) 



(1) 



where E is the neutrino energy and L is the flight length of neutrinos. Since the observed zenith 



angle distributions of e-like events agree with the predictions, 2-flavor f^ 



^-> v.. 



oscillation is 



considered to be dominant in the atmospheric neutrino oscillations. A x^ test is employed to 
evaluate the agreement of the fit to the observed data. The FC and PC samples are divided 
into 10 bins equally spaced between cos 6 = —1 and cos© = +1, where @ is zenith angle of 
particle direction. Furthermore, the FC sample is divided by the momentum, number of rings 
(single-ring and multi-ring) and particle types (e-like and /x-like). The upward stopping and 
upward through-going muon samples are divided into 5 and 10 bins equally spaced between 
cos© = — 1 and cos© = 0, respectively. In total 175 bins are used in the zenith angle analysis. 
A global scan was carried out on a (sin^ 29, Am^) grid including the unphysical region 
(sin^ 26 > 1). The best-fit for 2-flavor u^ <-5- Uj- oscillations was obtained at (sin^ 29 = 1.0, Am? = 
2.0 X 10-3 eV^) in which Xmin = 170.8/172 DOF. Figure El shows the contour plot of the 
allowed neutrino oscillation parameter regions. Three contours correspond to the 68%, 90% 
and 99% confidence levels (C.L.), respectively. The allowed oscillation parameter regions were 
sin^ 29 > 0.90 and 1.3 x lO^^ < Am^ < 3.0 x lO^^ eV^ at 90% C.L. 
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Figure 2: Allowed oscillation parameter regions for f^ <-» Vt oscillations. Three contours correspond to the 68 % 
(dotted line), 90% (solid line) and 99% (dashed line) C.L. allowed regions, respectively. 



2.3 L/E analysis 

Zenith angle dependent deficit of muon neutrinos has been interpreted as evidence for atmo- 
spheric neutrino oscillations l^^. Neutrino oscillations have been studied using various neutrino 
sources. However, the sinusoidal neutrino flavor transition probability predicted by neutrino 
oscillation has not been demonstrated yet. The analysis described herein used a selected sample 
of atmospheric neutrino events, those with good resolution in L/E, to search for the dip in 
oscillation probability expected when the argument of the second sine-squared term in Eq. ^ is 
7r/2. 

In the L/E analysis, 1489 live-days exposure of FC //-like and PC atmospheric neutrino data 
are used. Event selection and classification in L/E analysis are different from those in the zenith 
angle analysis. In order to increase the statistics of the data, especially of high energy muons, 
the fiducial volume for the FC sample is expanded from 22.5 kton to 26.4 kton. Estimated non- 
neutrino background events in the expanded fiducial volume is less than 0.1 % and is negligibly 
small. The PC events is subdivided into two categories : "CD stopping events" where the muon 
stops in the outer detector, and "OD through-going events" where the muon exit into the rock. 
The division is based on the amount of Cherenkov light detected in the OD. Since these two 
samples have different resolution in L/E, different cuts were applied for each sample, improving 
the overall efficiency. 

The neutrino energy is estimated from the total energy of charged particles observed in the 
ID. The energy deposited in the OD is estimated from the potential track length in the OD 



and is taken into account for PC events. The relationship between the neutrino energy and 
the observed energy is determined based on the Monte Carlo simulation. The flight length of 
neutrinos, which ranges from approximately 15 km to 13,000 km depending on the zenith angle, 
is estimated from the direction of the total momentum of the observed particles. 
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Figure 3: Contour plots of 70 % L/E resolution on the (cos O, E^) plane for (a) FC single-ring, (b) FC multi-ring, 

(c) PC OD stopping and (d) PC OD through-going samples. Three lines in the upper left figure show the survival 

probabilities of muon neutrinos predicted from neutrino oscillation with (sin'^ 26, Am^) — (1.00, 2.4 x 10~^ eV^). 

Full and half oscillation occur on the solid and dashed lines, respectively. 



The resolution of the reconstructed L/E is calculated at each point on the (cos 0, Ey) plane, 
where is the zenith angle. Figure El shows 70 % L/E resolution contours which are used as 
the selection criteria. The L/E resolution cut is set to be A{L/E) < 70% from the Monte 
Carlo simulation to maximize the sensitivities to distinguish neutrino oscillation from other hy- 
potheses. Three lines in Figure |31 (a) indicate survival probabilities of muon neutrinos predicted 
from neutrino oscillations. Statistics of high energy muon events is crucial to observe the first 
oscillation minimum in L/E. 

The left-hand plot in Figure |1] shows the number of events as a function of L/E for the data 
and Monte Carlo predictions, and the right-hand plot shows the data over non-oscillated Monte 
Carlo ratio with the best-fit expectation for 2-flavor f^ <-^ z^,- oscillations in which systematic 
errors are considered. A dip, which should correspond to the first oscillation minimum, was 
observed around L/E = 500 km/GeV. 

A fit to the observed L/E distribution was carried out assuming neutrino oscillations. In 
the analysis, the L/E distribution is divided into 43 bins from log{L/E) = 0.0 to 4.3 . The 
likelihood of the fit and the x^ are defined as : 
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where Nf^^ is the number of the observed events in the z-th bin and Nf^ is the number of 
predicted events, in which neutrino oscillation and systematic uncertainties are considered. 25 
systematic uncertainties are considered in the L/E analysis, which include uncertainty param- 
eters from the neutrino flux calculation, neutrino interaction models and detector performance. 
Among these, only 24 constrain the likelihood as the absolute normalization is allowed to be free. 
The second term in the likelihood definition represents the contributions from the systematic 
errors, where ctj is the estimated uncertainty in the parameter ej. 

A scan was carried out on a {siv?29, logA?n,^) grid, minimizing x^ by optimizing the sys- 
tematic error parameters. The Xmm value was 37.9/40 DOF obtained at (sin^20, /^m?) = 
(1.00, 2.4 X 10~^eV^). Including unphysical parameter region {s\v?26 > 1), the best-fit was 
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Figure 4: Left: Number of events as a function of the reconstructed L/E for the data (points) and the atmospheric 

neutrino Monte Carlo events (histogram). Right: Ratio of the data to the non-oscillated Monte Carlo events 

(points) with the best-fit expectation for 2-flavor v^ <-> Ut oscillations (solid line). Also shown are the best-fit 

expectation for neutrino decay (dashed line) and neutrino decoherence (dotted line). 
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Figure 5: 68, 90 and 99% C.L. allowed oscillation parameter regions for 2-flavor i/^ 

the L/E analysis. 



v-r oscillations obtained by 



obtained at (sin 26, Am,^) = (1.02, 2.4 x 10~^eV^), in which the minimum y^ was 0.12 lower 
than that in the physical region. Figure [S] shows the contour plot of the 68, 90 and 99% C.L. 
allowed oscillation parameter regions. The 90% C.L. allowed parameter region is obtained as 
1.9 X 10^3 eV^ < Am^ < 3.0 x 10"^ eV^ and sin^ 20 > 0.90. The location of the allowed region 
is consistent with that of the zenith angle analysis. ^^ 

The observed L/E distribution was also fit assuming neutrino decay'^^'^ and neutrino de- 
coherence!^ . The u^ survival probability for neutrino decay is expressed as P^iy^i — >■ Vfj) = 
[sin 9 + cos^ 0exp (— m/2T • L/E)\ where r is the lifetime of a neutrino mass state, and that 



for neutrino decoherence is expressed as P{y^ 



l — \ sin^ 26 [1 - exp {-^qL/E)] where 70 



is the decoherence parameter. The right-hand plot in Figure |1] includes the best-fit expectation 
for neutrino decay and decoherence. The \^in values were 49.1/40 DOF at (cos^0, m/T) = 
(0.33, 1.26 x 10-2 GeV/km) for neutrino decay and 52.4/40 DOF at (sin2 26', 70) = (1.00, 1.23 x 
10^2^ GeV/km) for neutrino decoherence. These Xmm values are 11.3 (3.4 standard deviations) 
and 14.5 (3.8 standard deviations) larger than that for neutrino oscillation. Alternative models 
that could explain the zenith-angle and energy dependent deficit of muon neutrinos are disfa- 
vored, since they do not predict the dip structure in the L/E distribution. 

The observed L/E distribution, especially the first dip, gives the first direct evidence that 



the neutrino flavor transition probability obeys the sinusoidal function as predicted by neutrino 
flavor oscillations. 

3 Solar neutrinos 

3.1 Introduction 

The observed solar neutrino flux in all the experiments has been significantly smaller than that 
was expected. From the recent solar neutrino data, especially together with Super-Kamiokande 
and SNO, the cause was identified to be neutrino oscillations'^. In order to determine the 
neutrino oscillation parameters, it it important to measure not only solar neutrino flux but also 
the energy spectrum and the time variations, which are independent of the uncertainties in the 
solar models. 

Solar neutrinos are observed in Super-Kamiokande with a large statistics owing to its large 
target volume. The detector can measure the energy spectrum of neutrino events precisely 
by the well calibrated performance^. Since the observation is working in real-time, the time 
dependence of solar neutrino flux, day/night or seasonal differences, can be also measured. In 
this paper, the results from 1496 days of solar neutrino data observed from May in 1996 to July 
in 2001 during the first phase of Super-Kamiokande are reported. 

3.2 Results of solar neutrino observations 

The solar neutrino signals observed in Super-Kamiokande are separated from background events 
by fitting the forward peak to the solar direction. The observed solar neutrino fiux in Super- 
Kamiokande, whose energy threshold is 5.0MeV, is : 

2.35±0.02(stat.)±0.08(sys.) [xlO'^ /cm'^ /sec]. (4) 

As compared with the result to standard solar model (BP20007^, the ratio is : 

^ = 0.465 ± 0.005{stat.n'olt{sys.), (5) 

The observed flux is significantly smaller than the prediction. The day-night fiux differences are 
also observed. A simple day/night asymmetry is obtained by dividing the data sample into day 
and night as follows : 

^DN = ,/'V/"^'w9 = -0-021 ± 0.020{stat.)1:'of,l{sys.), (6) 

[9 day + 9night)/^ 

which is consistent with zero. 

3.3 Solar neutrino oscillation analysis 

Maximum likelihood fit was carried out to determine the a llow ed neutrino oscillation parameters 
considering the time variations of the solar neutrino fiuxl^. The data sample is divided into 
21 energy bins from 5 MeV to 20 MeV. Two types of probability distribution functions are 
considered. p{cos 6aun, E) describes the angular shape expected for solar neutrino signals of 
energy E and tij(cos0sun) is the background shape in the energy bin i, where ^sun is the angle 
between the reconstructed recoil electron direction and the solar direction. The likelihood : 

jr = e-iT..B.+s) -Q Yl lBi.Ui{cosei^) + S—^-p{cosei^,E^)-z{a,t,)\ (7) 

1 = 1 K = l \ ^j ■? / 



=0.75 




Figure 6: Observed spectrum (top) and D/N asymmetry (bottom) . The predictions (solid lines) are for tan'^ 9 = 

0.55 and Am^ = 6.3 x lO'^^eV^ with 08^ = 0.96SSM and 4>hep = 3.6SSM. Each energy bin is fit independently to 

the rate (top) and the day/night asymmetry (bottom). The gray band is the ilcr range. 

is maximized with the respect to the signal S and the 21 backgrounds Bi. MCj is the number 
of events expected in the energy bin i. To take into account the time variations of the solar 
neutrino flux in the likelihood fit, the signal term is modified by z{a,t^), where t^ is the event 
time and a is an amplitude scaling factor. Solar zenith angle (day /night) variations and the 
additional seasonal variation are considered. For the best-fit LMA parameters, the amplitude 
scaling factor is obtained as a = 0.86 ± 0.77, which corresponds to the day/night asymmetry : 



Aj3N = -0.018 ± 0.016(stat.) 



-0.013 
-0.012 



(sys.) 



(8) 



where —2.1% is expected from the simple day/night comparison. The statistical uncertainty 
is reduced by 25 % with this likelihood analysis. However, the resulting day/night asymmetry 
is still consistent with zero. Figure IHl shows the resulting fitted rate (top) and the day/night 
asymmetry (bottom) for each energy bin with the expectations for the best-fit LMA oscillation 
parameters. 

To constrain neutrino oscillation parameters using the time variations of solar neutrino 
flux, the likelihood difference Alog>C = log>C(Q; = 1) — log C{a = 0) between the expected 
time variation and no time variation is computed. AlogjC is interpreted as a time- variation 
Ax^y = — 2Alog>C and is added to the spectrum x^- The total x^ is defined as : 
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where d, is the ratio of the data to the SSM prediction in the energy bin i and pi is the prediction 
in which neutrino oscillation is considered, at is energy bin-uncorrelated uncertainty, as, <7s 
and aji are the estimated uncertainties in the ^B neutrino spectrum, SK energy scale and SK 
energy resolution, respectively. The last term constraining the ^B flux to the SSM is optional. 
The left-hand plot in Figure [7| shows the allowed and excluded regions at 95% C.L. obtained 
from the SK solar neutrino data. The advantage of this method is that the number of energy 
bins can be increased because the analysis is independent of the binning in night time. 

Stronger constraints on neutrino oscillation parameters result from the combination of the 
SK measurements with other solar neutrino datal^ and the KamLAND results'^. The 95 % 
C.L. allowed region is shown in the right-hand plot in Figure [7| Only LMA solution remains. 
The allowed range of oscillation parameters are obtained to be Am^ = 7.21q'5 x lO^^eF^ and 



tan^ 



0.38 ±0.08. 
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Figure 7: Left: Excluded (SK spectrum and time variation; dark gray) and allowed (SK spectrum, rate and time 
variation; light gray) neutrino oscillation parameter regions at 95% C.L. Overlaid are the areas excluded by the 
day/night and seasonal variation (hatched regions). The plots at the top (right) show the x'^ difference as a 
function of tan^ 6 {Ani^). Right: Allowed region at 95% C.L from all solar data, KamLAND results, and the 
combination of all solar and KamLAND data. The plots at the top and right show the x^ differences. 



3.4 Summary 

Zenith angle analysis and L/E analysis were carried out using atmospheric neutrino data. A dip 
in the L/E distribution was observed for the first time, as predicted from the sinusoidal flavor 
transition probability of neutrino oscillation. The allowed neutrino oscillation parameter region 
was constrained to be 1.9 x lO'^ eV^ < Am^ < 3.0 x lO'^ eV^ and sin^ 29 > 0.90 at 90 % C.L., 
which is consistent with that of the zenith angle analysis. 

Super-Kamiokande has precisely measured the solar neutrino flux, recoil electron spectrum, 
and time variations of the flux. No significant time variation and energy distortion appeared 
from the 1496 live days of data. The time variation of neutrino flux was considered in neu- 
trino oscillation analysis by a maximum likelihood method. Combined with all solar neutrino 
experiments and KamLAND data, neutrino oscillation parameters are obtained to be Am? = 
7.2l[J:^ X lO-^eF^ and tan^ 6 = 0.38 ± 0.08. 
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